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O Abstract 

•i-H 

Population growth is in general constrained by food production, which 
4^ in turn depends on the access to water resources. At a country level, 

■ ^ some populations use more water than they control because of their 

ability to import food and the virtual water required for its production. 
^ Here, we investigate the dependence of demographic growth on avail- 

00 ^ble water resources for exporting and importing nations. By quanti- 

^ fying the carrying capacity of nations based on calculations of the vir- 

^ tual water available through the food trade network, we point to the 

* existence of a global water unbalance. We suggest that current export 

Q rates will not be maintained and consequently we question the long-run 

CO sustainability of the food trade system as a whole. Water rich regions 

are likely to soon reduce the amount of virtual water they export, thus 
^ ^ leaving import-dependent regions without enough water to sustain their 

^ populations. We also investigate the potential impact of possible scenar- 

^ ios that might mitigate these effects through (1) cooperative interactions 

among nations whereby water rich countries maintain a tiny fraction of 
their food production available for export; (2) changes in consumption 
patterns; and (3) a positive feedback between demographic growth and 
technological innovations. We find that these strategies may indeed re- 
duce the vulnerability of water-controlled societies. 
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1 Introduction 



Most of the water we use is to produce the food we eat. With the world's pop- 
ulation that has doubled every 40 years there is a growing concern that water 
limitations will soon impede humanity to meet its food requirements (1-7). The 
urgency to deal with this alarming situation by developing durable socio-political 
and economic strategies that promote a sustainable use of the environment and its 
natural resources, was at the core of the recent Rio-i-20 Earth Summit organized by 
United Nations (7-8). In recent years a number of studies have combined projec- 
tions of population growth with predictions of water availability and agricultural 
productivity under a variety of climate change and land use scenarios. These pre- 
dictions have been used to assess whether mankind will run out of water in the 
next few decades, and to investigate possible strategies to deal with the global 
water and food crisis (9-15). At the regional scale, there exist several areas of the 
world where the demand has already exceeded the supply of renewable freshwa- 
ter resources (16-19). How can this negative water budget be sustained? Mainly 
by importing food. The import of food commodities is associated with a virtual 
transfer of freshwater resources from production to consumption areas (20). Vir- 
tual water trade allows some populations to exceed the limits imposed by their 
local water resources (14, 21-22). By sustaining demographic growth above the 
regional carrying capacity, virtual water trade has mitigated the effects of drought 
and famine in many regions of world (14,23). Thus, the redistribution of virtual 
water (VW) resources often appears as a remedy to regional water crises (22-23). 
Even though, presently, most exporting countries can afford to sustain VW ex- 
ports, their demographic growth might soon limit the amount of VW resources 
they can place on the global market. At some point these societies will have to de- 
cide whether they want to sustain the existing export rates, or prefer to reduce the 
exports in order to meet their own food demand. What are the global implications 
of these two scenarios? To answer to this question, we need to relate demographic 
growth both to the available freshwater resources and to global patterns of VW 
trade. We first investigate current trends of demographic growth and then develop 
model-based predictions of how population is expected to change as water rich 
countries start reducing their exports. 
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2 The Carrying Capacities of Nations 



We classify all countries around the world into five groups, depending on their 
supply and demand of VW, and on the resulting balance or unbalance between 
available and consumed water resources (see Figure 1). We observe that the 
water-rich regions are in North and South America, Australia, and the former 
Soviet Union (or "eastern block"). These regions are known for being major VW 
exporters (26,27). Virtual water dependent regions (i.e., regions that need VW 
imports to meet their demand) are mainly in Europe, Mexico, and the western 
side of South America. Despite VW trade, large parts of Africa and Asia remain 
affected by water stress. Because food production is the major form of freshwater 
consumption by human societies (24-25), we calculate the "carrying capacity" of 
a nation, i.e. its maximum sustainable population, based on the water resources 
currently available for agriculture and livestock. The carrying capacity, though 
difficult to quantify, is a key notion for characterizing the relation existing be- 
tween demographic dynamics and their possible resource limitation (7-9). In here 
we will show that for almost one third of all world nations (i.e. water rich and 
VW dependent countries), the carrying capacity depends on their food availability, 
which, in turn, depends on the available water resources. Therefore, a quantitative 
estimate of the average local carrying capacity Kj^^ of country i is obtained di- 
viding the total water currently available for food production in that country (i.e., 
the current footprint of crops, grazing and livestocks (20,21,26)), by the volume 
of water, W^, used to produce the food consumed on average by one individual in 
that nation (16,19,20). Kl^^ is the maximum population sustainable with the avail- 
able local freshwater resources of country i and is expressed in terms of number 
of individuals. The virtual carrying capacity Ky is the maximum sustainable pop- 
ulation of country j when VW imports and exports are accounted for. Thus, the 
virtual carrying capacity depends on the structure of the global VW trade network 
(26,27). We also may incorporate in our modeling approach the impact on K of 
changes in consumption patterns, crops expansion, and increase in the efficiency 
of agricultural production (26). 

2.1 Relation between demographic growth and water availa- 
bility 

We analyze and model the relation between demographic growth and the water 
balance of water rich and VW dependent countries (for a total of N=52 nations) 
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using a stochastic logistic model. The population dynamics, of nation i are 
expressed as (26) 



where i = 1, 2, ...A^. is treated as a random variable to account for stochas- 
tic fluctuations in available resources and for uncertainties in their estimate. Thus, 

is expressed as = + ^ is a white Gaussian noise with mean (<^) = 
and covariance {^{t)^{s)) = (j]^i5{t — s) (26). (•) denotes the average with re- 
spect to the stochastic fluctuations, while a\-i is the intensity of the fluctuations, 
and 5(.) the Dirac's delta function. The model is based on two key parameters: 
and K\ The mean population growth coefficient 6i\ is determined for each nation 
using demographic data for the period of record 1970-2010 (26-28). To assess 
whether the demographic growth of a nation is driven by either local or virtual 
water availability, we consider both the case of population dynamics controlled 
by local water resources (i.e., with average "carrying capacity", = ^/oc) 
the effect of VW trade (i.e., K' = K'y). 

Figures 2-3 summarize the results of this analysis. We stress that the effective 
carrying capacities are obtained from direct evaluation of the actual virtual water 
availability and water footprint calculations (26) and not from fit of demographic 
data. The average demographic growth of water-rich countries is well described 
by the logistic model (Eq. (1)) with = Kl^^, which indicates that the growth 
of these populations is driven by local water availability and is not limited by 
the important export of VW sustained by these nations. Interestingly, several na- 
tions seem to have just reached a critical phase in their demographic growth: the 
point of separation between the local water and the VW regime (see (26)). On 
the other hand, demographic growth of VW trade dependent nations follows the 
logistic model with = Ky\ thus, VW inputs may sustain the demographic 
growth of populations that increasingly rely on external VW resources (26). In 
these countries population growth would not be sustainable without the import of 
VW from other regions. The fact that demographic growth in water-rich countries 
is independent of their VW exports, while water-poor nations increasingly rely 
on VW imports highlights a situation that, in the long run, will be unbalanced 
and unsustainable. These results indicate that both water-rich and trade dependent 
populations are growing relying on the same pool of resources. At some point 
it will happen that, to meet their own demand, water-rich countries will have to 
reduce their exports, thereby causing the emergence of water limitations in trade 
dependent countries. Unless new freshwater resources become available or invest- 




(1) 
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ments for a more water-efficient agriculture are made, these populations will have 
to decrease. 

3 Sustainability of Future Scenarios 

To investigate these effects we use the logistic model (Eq. (1)) coupled with the 
dynamics of food trade through the VW network (Figure 3a). We first build a bi- 
partite network with nodes that interconnects m water rich nations (blue nodes 
r = 1, 2, ..m) to VW trade dependent countries (red nodes p = 1, N — m). 
We investigate two scenarios. In the first one,the population in blue nodes grows 
following Eq. ([T]) with = {1 — /3)Ki^^, where /5 is the maximum fraction 
of Ki^^ that water-rich countries are willing to place in the long term (t oc) 
in the global VW trade market. At any time t, a water rich country, r, with 
population x^{t) shares with trade dependent nations a volume of water equiv- 
alent to Ki^^ — we assume that this water is equally partitioned among 
all the dr trade dependent countries connected to that blue node (i.e., dr is the 
degree of node r). Therefore the population of a red node, p, grows accord- 
ing to the Eq. ^ with virtual carrying capacity = Kf^^ + XlJ^i ^rp[Kl^c ~ 
x'^{t))Wl/ {W^dr), where a^^ is the adjacency matrix describing the virtual water 
trade network (a^p=l if r and p are connected; a^p=0, otherwise). The case /3 = 
corresponds to the purely competitive case, where in the long run K[^^ 
and no water is exported because all resources are used to support the local pop- 
ulation of the blue node r. The second scenario takes into account plausible 
crop expansions, increases in agricultural production efficiency, and changes in 
diet and consumption rates (19,23,26), while the cooperation regime is turned off 
(/3 = 0). The relevant dynamics are expressed by the same equations as in the pre- 
vious case, but critically endowed with time dependent local carrying capacities, 
Kioc{t), and per capita virtual water consumption, Wc{t). In particular, we impose 
a growth rate in the local resources (i.e., in Kioc{t)) to account for possible crop- 
land expansions and crop yield enhancement in existing croplands, and a decrease 
in per capita virtual water consumption (i.e., in Wc{t)) due to new practices that 
reduce water loss (19) and changes in consumption patterns (29). 
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4 Discussion and Conclusions 



To evaluate the impact of the network structure, we investigate the system's dy- 
namics in the two above scenarios using four types of networks with different 
structural properties but the same number of nodes and average number of links. 
For the first scenario, we find that in the purely competitive case (f3 = 0) the pop- 
ulation of trade dependent countries increases in the first 25 years and decreases 
in the subsequent years as a result of water limitations arising from the exclusion 
from the access to water resources controlled and claimed by water-rich countries 
(26,30). This decline continues in the following decades as a results of demo- 
graphic growth in water-rich populations (Figure 3b). This behavior appears to be 
very robust and independent of the topological properties of the underlying net- 
work. If trade dependent countries use more VW resources than they control, their 
population will strongly decrease once the water rich countries start reclaiming 
all the water resources they have access to, no matter how water is redistributed 
among the trade dependent virtual water countries. The results are robust and do 
not qualitatively change with different strengths of the fluctuations of K in the 
range < cr^ < 0.20^. Interestingly, this decrease in trade-dependent pop- 
ulation is reduced if a cooperative regime (i.e., ^5 > 0) is considered. If water 
rich countries keep a fraction, /? > 0, of their water resources in the VW global 
market, VW dependent countries can sustain a larger population, which increases 
as a function of f3 (Figure 3b). The overall effect of a cooperative regime is a 
long-term increase in the total global population and thus a more sustainable de- 
mographic growth. In this cooperative regime, the network topology affects the 
coupled VW-demographic dynamics (26). We finally note that the sensitivity of 
the demographic dynamics on f3 is very strong, and just a small departure of /? 
from zero, may lead to substantial reductions in the decrease of trade dependent 
populations (Figure 3b and (26)). Also the second strategy, based on the enhance- 
ment of productivity efficiency and a decrease in per capita global consumption, 
results in a remarkable relief for trade dependent countries, whose populations are 
subjected to less pronounced declines (Figure 3c). We also find that the increase in 
food (and virtual water) availability resulting only through changes in consump- 
tion patterns and greater equity in per capita consumption would not be sufficient 
anyhow to meet the increasing demand of a growing human population. In fact, 
in the long run, if the growth rate of Kioc{t) tends to zero, the trade dependent 
populations will peak and then inevitably decline (26). 

Despite the presence of a number of other environmental, cultural, and health 
related factors not included in this study, this analysis points out how VW trade 
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is only a temporary solution to a local-to-regional unbalance between populations 
and food production. The existence of this unbalanced condition might be miti- 
gated if a cooperative regime among water rich and VW dependent nations will 
continue to exist even once the excess of VW in the exporting countries is strongly 
reduced by their demographic growth (/? > 0). We finally show that strategies 
aiming at an increase in productivity efficiency through agricultural practices that 
enhance crop yields while reducing water losses (e.g., water harvesting, water 
conservation, genetically modified crops) improve the sustainability of trade de- 
pendent societies with respect to a decrease in export rates from water rich coun- 
tries. 
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^ water rich nations 
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^ VW dependent 
O barely self-sufficient 
O water scarce 
O not significant data 



Figure 1: Map of the World's nations classified on the basis of their dependency 
on local and virtual water resources using data for the 1996-2005 period (20,21). 
Countries are water rich when their mean population, x, is less than CS^/^c; 
virtual water dependent if Ky > x > Kioc', barely self sufficient if Kioc ^ 
water scarce if x > Ky > Kioc Countries for which data do not appear to be 
consistent are shown in gray. A separate analysis has been carried out for the 
countries from the influence zone of the former Soviet Union (or the "eastern 
block") because in the past two decades their demographic dynamics have been 
affected by major political changes not related to water (26). 
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Figure 2: Comparison between the logistic model given by Eq. (1) and demo- 
graphic data (red dots) from 1970 to 201 1 (28) for three water-rich (blue) and four 
virtual water dependent countries (red). These nations are connected through the 
VW trade web. The green dashed lines represent the numerical simulations of 
Eq. (1) with carrying capacity = K^^^, while the blue lines corresponds to 
for virtual carrying capacity = Ky. Figures of the results for the complete 
analysis can be found in (26). The evolution of water rich countries population is 
predicted by Eq. (1) with = Kj^^, whereas for water dependent countries the 
appropriate evolution is given by Eq. (1) with = Ky. 
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Figure 3: a) Simulated coupled dynamical network. Blue nodes represent water 
rich nations, while red nodes are the VW trade dependent countries. We investi- 
gate the logistic model of coupled population dynamics on four different kind of 
networks: Random graph (26) (blue lines), small world network (26) (red lines), 
scale free network (26) (gray lines) and (in green) the graph with topological 
properties similar to those observed in the real global VW trade network (27) . b) 
Total population Xy^d of all VW dependent nations as a function of time in the 
pure competitive case (/3 = 0, solid lines), and in two different cooperative regime 
(/3 = 0.04-dot lines, and /3 = 0.08 dashed lines), c) Demographic dynamics of 
VW dependent nations according to scenario 2 (solid lines), which accounts for 
an increase in productivity efficiency, crop expansion and decrease in consump- 
tion rates in the absence of cooperation (/3 = 0). The color scheme is the same as 
in panel b. The dashed line corresponds to scenario 1 (panel b) with the existing 
virtual water trade network (27) and with no cooperation, /3 = 0. 
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